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Abstract 

In light of recent analyses that have shown that nontrivial hidden-sector dynamics in 
models of supersymmetry breaking can lead to a significant impact on the predicted 
low-energy supersymmetric spectrum, we extend these studies to consider hidden-sector 
effects in extensions of the MSSM to include a seesaw model for neutrino masses. A 
dynamical hidden sector in an interval of mass scales below the seesaw scale would yield 
renormalization-group running involving both the anomalous dimension from the hidden 
sector and the seesaw-extended MSSM renormalization group equations (RGEs). These 
effects interfere in general, altering the generational mixing of the sleptons, and allowing 
for a substantial change to the expected level of charged-lepton flavour violation in seesaw- 
extended MSSM models. These results provide further support for recent theoretical 
observations that knowledge of the hidden sector is required in order to make concrete low- 
energy predictions, if the hidden sector is strongly coupled. In particular, hidden-sector 
dynamics may impact our ability to reconstruct the supersymmetric seesaw parameters 
from low-energy observations. 



1 Introduction 



Weak-scale supersymmetry in the guise of the minimal supersymmetric extension of the 
Standard Model (MSSM) [TJ [2] provides an elegant solution to the gauge hierarchy prob- 
lem, satisfies electroweak precision constraints, and predicts weakly-interacting dark mat- 
ter (see, for example, [3]). However, the MSSM requires additional physics responsible 
for breaking supersymmetry itself (for reviews see [H EJ [6], [7]). This additional physics 
is expected to include a hidden sector that breaks supersymmetry spontaneously, and 
a messenger sector, which communicates the symmetry breakdown to the visible sector 
fields of the MSSM. The MSSM predicts gauge-coupling-constant unification [HJ [9j [10] 
at My ~ 2 x 10 16 GeV if there is a desert between the weak scale and the unification 
scale. The success of this prediction encourages the hope that one may gain insights into 
high-scale physics through the renormalization group equations (RGEs) of the MSSM. 

The observation of neutrino oscillations [TT1 [T^ 1131 [131 IISl IISl IIZl IIHl IISl I2DI Ell [231 [331 
EMI [23| I21H |2~T] suggestions the existence of new high-scale physics responsible for the small 
neutrino masses inferred from experiment. A popular framework for generating small 
neutrino masses invokes heavy Majorana gauge-singlet fermions that, once integrated 
out of the theory, yield the dimension-five operator LLHH with a coefficient suppressed 
by a factor of the heavy gauge-singlet mass scale. This framework - called the seesaw 
mechanism (for a review, see |28j) - in its simplest incarnation (the Type-I seesaw), 
demands that the Majorana mass scale appears below the unification scale, in order to 
satisfy the constraints provided by neutrino oscillations, unitary, and triviality [2"5l I5UI 
EHE2], and typical values of the Majorana scale are near Mr ~ 10 14 GeV. 

Since the seesaw mechanism violates lepton number by two units, its natural super- 
symmetric extension is consistent with R-parity conservation, preserving the stability of 
the lightest supersymmetric particle and hence its candidacy for cold dark matter. The 
supersymmetric seesaw can be therefore be incorporated as a simple extension of the 
MSSM. The appearance of the Majorana scale in the desert below the unification scale 
implies an interval in which the RGEs of the MSSM are augmented by the interactions 
contributing to the seesaw mechanism. The change in the renormalization-group flow of 
the MSSM induced by the heavy gauge-singlet Majorana fermions does not affect the 
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RGE evolution of the gauge couplings at leading order, but does, in general, generate 
off-diagonal mixing in the slepton sector that can lead to a observable predictions for 
charged-lepton flavour-violating decays. This prospect offers further insight into high- 
scale physics from low-energy observations. 

While it has been long understood that different means of supersymmetry breaking 
lead to different sparticle spectra and different low-energy phenomenologies [H E31 EU 
[351 [361 EZl EH [391 SOI SH S21 S3l [H], only recently has it been realized that hidden- 
sector dynamics may play an important role in low-energy predictions [45l 1461 1471 [48]. 
Specifically, it has been shown that if the hidden sector contains strong self-interactions, 
hidden-sector renormalization effects can influence significantly the renormalization-group 
running of the MSSM scalar sector. These hidden-sector renormalization effects may 
modify in an observable way the simple mass relations expected naively from scalar-mass 
unification at the unification scale. 

Observable-sector effects from the hidden sector result from quantum corrections that 
correct the non-renormalizable operators introduced at the messenger scale required to 
communicate supersymmetry breaking to the observable sector. Since the hidden-sector 
scale sits in the desert below the unification scale (e.g., Madden ~ 10 12 GeV in typical 
gravity- mediated models), hidden-sector renormalization effects will be present along- 
side the seesaw mechanism. Furthermore, since the hidden-sector renormalization af- 
fects directly the diagonal scalar-mass terms, the effect will alter the ratio between the 
radiatively-generated seesaw off-diagonal slepton masses and the diagonal terms. In gen- 
eral, models with moderately strongly-coupled hidden sectors will therefore have an im- 
pact on the expected amount of charged-lepton flavour violation. 

This paper is a continuation of our previous work [49J in which we explored the 
possible measurability of the diagonal scalar-mass effects. Here we explore the effects 
of hidden-sector renormalization on the seesaw extension of the MSSM with gravity- 
mediated supersymmetry breaking. We examine the consequent impact of hidden-sector 
effects on charged-lepton flavour violation, in particular on the induced rate for fi — > 
in the seesaw extension of the MSSM, using model-independent anomalous- dimension 
parametrizations of the hidden sector's renormalization group evolution. In this fashion, 
we examine a wide range of model classes for the hidden sector and their effect on the 
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predictions for charged-lepton flavour violation in models with a supersymmetric see- 
saw. We also discuss the impact of these effects on the possibility of reconstructing the 
supersymmetric seesaw parameters from weak-scale observations. 

2 The MSSM Seesaw 

We consider the class of gravity-mediated supersymmetry-breaking models that lead to 
the constrained MSSM (CMSSM) at the unification scale, with universal flavour-diagonal 
soft masses, universal gaugino masses, and tri-linear A-terms proportional to the superpo- 
tential Yukawa couplings and the universal scalar mass. Charged-lepton flavour violation 
in supersymmetric seesaw models arises from renormalization-group running of the see- 
saw sector between the unification scale and the Majorana mass scale [50| [51]. The 
seesaw sector thereby induces radiatively off-diagonal slepton mass terms that contribute 
to flavour-violating decays. 

There is an elegant parametrization [52] for encoding the seesaw parameters. In 
the basis where the heavy neutrino singlets, the charged-lepton Yukawa matrix, and 
gauge interactions all appear flavour- diagonal, and where Y„ denotes the seesaw Yukawa 
couplings, the product YlY u can be written as, 

YtY u = U PM m^rtMR^Ul um , (1) 
where k contains the light neutrino masses inferred from low-energy experiments, 

My 

M u = diag(m ui , m^,™^) , (2) 

and M denotes the diagonal Majorana singlet mass matrix, diag (Mi, M2, Ms). The 
matrix £/pmns labels the neutrino mixing matrix inferred from the neutrino oscillation 
data, and the orthogonal matrix R contains, in principle, three additional complex mixing 
parameters originating at the Majorana mass scale. Eq.flTJ provides a general description 
of the seesaw mechanism and gives a useful parametrization for examining lepton-flavour 
violation (LFV) in seesaw models. 

In order to determine the level of LFV in a given model, the full MSSM RGEs must 
be integrated from the unification scale to the weak scale, integrating each Majorana 
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gauge-singlet neutrino out successively at its appropriate mass scale. In the following 
sections we examine two interesting limiting cases of the seesaw: strongly hierarchical 
Majorana gauge-singlet neutrinos and degenerate Majorana gauge-singlet neutrinos, both 
with hierarchical light neutrinos, and each with the Majorana scale set at Mr = 10 14 GeV. 
This restriction reduces the number of free parameters in the orthogonal matrix R. 

3 Hidden-Sector Renormalization Effects on LFV at 
Leading-Logarithmic Order 

As a demonstration of the basic idea, we consider the toy self-interacting hidden-sector 
presented in [46j|45], which contains the cubic superpotential term 



This simple superpotential cannot by itself break supersymmetry, and hence the hidden 
sector must contain additional interactions responsible for generating an F or D-term 
vacuum expectation value (VEV) in any realistic model. For the purposes of examining 
the effects of hidden-sector renormalization, in this section we will suppose that eq.([3]) 
appears as the dominant self-interaction term in the hidden-sector superpotential, and 
that it provides the dominant hidden-sector contribution to the anomalous dimension of 
the operator mediating supersymmetry-breaking scalar masses in the observable sector. 
We assume that the hidden-sector chiral superfields couple to the visible sector fields 
through the non-renormalizable operators 



where $j and W n denote the MSSM chiral superfields and gauge fields respectively, and 
Mpi denotes the scale of gravitational mediation - the reduced Planck mass. Once super- 
symmetry breaks in the hidden sector, these terms yield the usual soft supersymmetry- 
breaking masses in the MSSM. In particular, fcj generates the soft scalar mass terms, 
whilst u determines the gaugino masses. In holomorphic renormalization schemes, non- 
renormalization theorems protect to at all scales. Finally, the usual MSSM gauge and 
Yukawa interactions generate the MSSM RGEs. 



W h = -X 3 . 



(3) 
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Figure 1: One-loop supergraphs contributing to the scalar-mass renormalization. The first 
diagram represents the quantum corrections arising from the hidden sector, whilst the last 
two diagrams represent the usual visible- sector loop diagrams. 

The coefficient ki in eq.fll]), which we take to be flavour- diagonal at the unification 
scale, is renormalized by the two separate contributions given in Fig. [TJ the usual visible- 
sector interactions involving vector and chiral superfields of the MSSM, and the hidden- 
sector self-interactions in eq.(j3]). As a result, the RGE for ki in the presence of the hidden 
cubic superpotential reads 

dh 2A*A, , . , . 

-ki + gauge and visible sector contributions. (5) 



dt 16tt 2 

The hidden-sector renormalization implies that the scalar-mass RGEs are augmented 
between the hidden-sector and messenger scales. In the hidden-sector theory under con- 
sideration, the RGE for an MSSM scalar sparticle becomes 

dm s 2 dm s 2 2A*A 2 

^r^^r + i6^ ms - (6) 

In general, the additional terms augmenting the usual MSSM RGEs arising from the 
hidden sector will be more complicated than the prescription given in eq.flE])- To give a 
full description of the hidden-sector effect on the MSSM RGEs, a complete model of the 
hidden sector would be required [UJ 051 0SJ • We stress that in this section we are simply 
considering a toy example that illustrates the effects on the level of slepton mixing, and 
hence of charged LFV, in supersymmetric models. We return to the problem of model 
dependence in following sections. 

While eq.(E]) demonstrates that we may expect shifts in the mass eigenvalues in the 
scalar sector of the MSSM, the hidden sector can also shift indirectly the off-diagonal 
elements. Given that the off-diagonal scalar mass-matrix elements determine the level 



of flavour- changing neutral currents in the MSSM, it is worth examining the size of this 
indirect effect. In particular, given that the seesaw sector sits near the hidden-sector scale 
in the model class we consider, the effect on the off-diagonal structure arising from the 
hidden sector will impact the amount of predicted charged LFV. To get an idea of the 
underlying physics of the effect, we proceed with our toy analytical example with scalar 
masses-squared running augmented by hidden-sector effects, but for analytic simplicity 
and clarity we ignore for the moment the running of the hidden-sector coupling A itself, 
which is 0(A 3 ). 

In the MSSM with Majorana gauge-singlet neutrinos, the RGE of the left-handed 
slepton mass-squared reads (before including the hidden-sector contribution): 
<im L 2 



= m L ^Y e + Y e T Y e m L ^ + m^Yl Y u + Y\ Y„nu/ 



+2 (Ytm e 2 Y e + m 2 H Y\Y e + A^A e ) + 2 (Y^rn^Y, + m^Y+Y,, + A l A -) 



-(2g 2 1 \M 1 \ 2 + Qg 2 2 \M 2 \ 2 ) I 3 , 

where t denotes the logarithm of the running scale, and Y v denotes the seesaw Yukawa 
couplings. In order to examine the leading-logarithmic behaviour for the prediction of 
LFV in the presence of hidden-sector effects in this model class, we consider the following 
simplified 2-by-2 model case: 

r I 
I 



d 
dl 



2/n 2/12 

2/21 2/22 



A 



2/ii 2/12 

2/21 2/22 



+ 3 



n n n X2 
n 2 i n 22 



2/ii 2/12 

2/21 2/22 



(?) 



In eq.flEJ), we model the hidden-sector effect by the term proportional to A. In this case 
the matrix 

2/ii 2/12 

2/21 2/22 

is analogous to the slepton mass-squared matrix, and 



(9) 



(10) 



nn n X2 
n 2 x n 22 

is analogous to the matrix product Y\Y V . We place a factor of three in eq.flEl) to match 
the factor three that appears in the leading-logarithmic analysis of eq.(JZJ). In our ap- 
proximation we ignore the possible presence of A-terms. The initial conditions at the 
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unification scale read: 
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2/12 






ml 





2/21 


2/22 










ml _ 



(11) 



and we assume that yu,y 22 2/12,2/21 throughout. Note that the element y 2 \ controls 
the level for the branching ratio // — >• e7- Using the above matrix differential equation, 
we obtain a coupled set of differential equations, namely, 

dy 2 i 3 



dt 
dyu 



(^222/21 + ^212/11) 
((A 2 + 3n u )yn + 3ni 2 2/2i) • 



di 8tt 2 

In the limit that y n ^> y 2 \, the approximate solution for yu becomes 



(12) 
(13) 

(14) 



and using this solution, we obtain 

dy 21 3 



dt 



2 v ^-222/21 + n 2 im e K 



8tt 

3^21^0 (A 2 +3nn)/87r 2 t 
87T 2 



(15) 



where the approximation again makes use of ^> y 2 \ ■ Accordingly, the approximate 
solution for y 2 \ reads, 



y 21 = ( e (\ 2 +3n n )/SnH _ j 



(16) 



A 2 + 3n n 

Integrating the renormalization group flow to the seesaw scale yields t = — ln(Mx/MR 
which allows us to write 



_ 3n 21 m / (A 2 +3rtll)/87r 2 i n (M x /M R ) 
1/21 A 2 + 3n u V 

Two limiting cases of eq. ffTTl) are worth exploring: 



1 . 



(17) 



A 2 1 In this limit, we recover the usual leading- logarithmic result for seesaw- 
induced off-diagonal slepton mass terms, namely 



2/21 



3n 21 ml , ( M x 



8tt 2 



In 



Mi 



3(YtY,) 21 m 2 h)/ ■ 



where we have made the identification n 2 \ — > (YjYj,) in the last line above. 
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(18) 



A 2 1ii(Mx/Mr) > 87r 2 In this limit, assuming A ^> n u , we now have 

-A 2 /8tt 2 









\m r 


3 (yjy„; 


21 m 


A 2 



Mi 



R 



-A 2 /8tt 2 

- 1 I , (19) 



where again we made the identification n 2 \ — > (YjY^) . We see that the j/ 2 i 
depends on 1/A 2 , which can lead to a suppression in the mixing for large A. 

The branching ratio for fi — > e7 depends at leading order on the off-diagonal slepton 
mass matrix, mL 2 i, via 

3 I 2 I 2 

BR(/i -> e 7 ) « -^ |mL211 tan 2 /?, (20) 

where mj denotes a typical sparticle mass. The leading-logarithmic approximation for 
the off-diagonal slepton mass matrix elements yields 

ml os , / 



m 4 



where m 2 , denotes the universal scalar mass, a® labels the constant of proportionality in 
the A-terms, and Mr denotes the Majorana mass scale. The branching ratio for fi — > e 7 
therefore becomes 



BR(/i e 7 ) 



a 3 



G 2 F m 8 s 



■(3 + a ) In 1 



(YtY,) 2 irtan 2 /3. (22) 



8tt 2V ' \M R/ 

The current experimental bound reads BR(/i — > e 7 ) < 1.2 x 10 -11 , and the MEG exper- 
iment at PSI expects to attain a sensitivity at a level of BR(/i — > e 7 ) ~ 5 x 10~ 14 [53] . 
In the next section we calculate the branching ratio for /1 — > e 7 using the full one-loop 
expression arising from Fig. [2j and we run the full one-loop RGEs for the MSSM, in 
classes of parametrizations of the hidden sector. 

Using eq. ff2TT) . the branching ratio for fi — > e 7 depends on y 2 i as 

BR(/i -> e 7 ) ~ i||L tan 2 /3. (23) 

We see from Fig. [3] that the suppression of the branching ratio as a function of the 
fixed hidden-sector coupling becomes larger than a factor of two when A > 5. We 
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Figure 2: One-loop diagrams contributing to jj — > c~) in the MSSM. The symbols fb denote 
charged sleptons, denote sneutrinos associated with the left-handed slepton doublet, 
denote the charginos, and x° denote the neutralinos. 
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Figure 3: The branching ratio for [l — *> e7 calculated incorporating the hidden sector in 
leading-logarithmic order as a function of X, relative to the value without hidden-sector 
effects. 

should emphasize that in the above analysis we did not consider the shift in the scalar 
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spectrum relative to the MSSM from hidden-sector dynamics, nor did we include the 
renormalization-group running of A. Since the branching ratio depends inversely on spar- 
ticle masses to the eighth power (see eq. fl23l) ). the changes in the sparticle spectra will 
have a non-negligible effect on the predicted rate. A competition will emerge between 
the suppression expected from the above analysis and the shifts in the sparticle spectra. 
We anticipate that in some regions where the hidden sector lowers the scalar spectrum 
sufficiently relative to the unaltered MSSM, the light scalars will dominate over the sup- 
pression factor, and yield an enhanced rate for fi — > cy. We explore these details in the 
next section. 

Our analysis thus far assumed A was constant over the range of integration, in order 
to allow simple analytic treatment. In reality, the hidden-sector coupling A runs with its 
own (3 function: [3\ = (3/(327r 2 ))A 3 . However, since the superpotential of eq.® by itself 
does not break supersymmetry, we must have additional hidden-sector interactions, so 
this model is just a toy, and not to be taken literally in detail. The point of the above 
discussion was simply to illustrate how a hidden sector self-interaction can influence the 
level of predicted charged LFV. In the following section, we examine the impact of the 
hidden sector in a model-independent fashion. Instead of hypothesizing a particular form 
of the coupling in the hidden sector, we parametrize the anomalous dimension itself. This 
parametrization allows us to examine both IR- and UV-free theories with ease. 



4 LFV with general Hidden-Sector Effects 

We saw from eq.© that in the cubic hidden superpotential theory we had 

dh 2A*A, . , . 

= -ki + visible sector contributions. (24) 

dt 16vr 2 v ; 

In order to examine general theories of the hidden sector, we require a parametrization 

that does not depend on the particulars of the hidden-sector implementation. For a 

general hidden sector we can define 

dk- 

= j(t)ki + visible sector contributions, (25) 

dt 

where j(t) denotes the anomalous dimension contributed by the hidden sector. Instead 
of using 7(t) from some specific theory as we did in the previous section, following our 
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previous paper [49J, we consider the parametrization 

7(*) = j-r^ y (26) 

where fe 7 and a 7 are theory-dependent factors, and s denotes the logarithm of the running 
scale. Depending on the location of the pole in 7(t), we can examine both IR- and UV- 
free theories. We consider the following four cases for the anomalous dimension ~f(t): 
6 7 = ±2/3, a 7 = log(5 x 10 n GeV//i), log(5 x 10 18 GeV//i). We depict these cases in 
Fig.! 




Renormalization scale (GeV) Renormalization scale (GeV) 




Renormalization scale (GeV) Renormalization scale (GeV) 



Figure 4: Four panels demonstrating the possible behaviours of the anomalous dimension 
outlined in the text. Top left: b 1 = 2/3, a 7 = log (5 x 10 11 GeV/fi), top right: 6 7 = —2/3, 
a = log(5 x 10 11 GeV/fi), bottom left: 6 7 = -2/3, a = log(5 x 10 18 GeV//i), bottom 
right: 6 7 = 2/3, a = log (5 x 10 18 GeV//i). In each panel, the heavy solid line at 2 x 10 18 
GeV denotes the reduced Planck mass, the light solid line at 2 x 10 16 GeV denotes the 
unification scale, the heavy dashed line at 10 14 GeV represents the Majorana mass scale, 
and the light dashed-dotted line at 10 12 GeV represents the hidden-sector mass scale. 

By choosing the poles in either of two locations - just beyond the reduced Planck 
mass, or just below the hidden-sector scale (10 12 GeV) - we arrange that the hidden sector 
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remains perturbative over the integration range and up to the reduced Planck mass itself. 
Once we reach the hidden-sector mass scale, we integrate out the hidden-sector physics, 
returning to the usual MSSM RGEs. As we see in Fig. HI by placing the pole above the 
reduced Planck mass in the IR-free case, the magnitude of the anomalous dimension in the 
interval between the unification and the Majorana mass scales is smaller than in our UV- 
free case. We make this choice in order to ensure that the Landau pole in the IR-free case 
does not appear below the reduced Planck mass. In the following section, we will apply 
each of these cases to the MSSM with the Majorana scale set at 10 14 GeV and with both 
hierarchical and degenerate Majorana gauge-singlet neutrinos, under the assumption of a 
normal hierarchy for the light neutrinos of the Standard Model. To determine the effect 
of the hidden sector on the predicted rate for \i — > ej, we run numerically the seesaw- 
extended MSSM RGEs including the hidden sector from the unification scale to the weak 
scale. We integrate out the Majorana gauge-singlet neutrinos at their associated scale 
and we also integrate out the hidden-sector effect at the hidden sector scale of 10 12 GeV. 

4.1 Hierarchical Gauge-Singlet Neutrinos 

We recall from section [2] that the branching ratio BR(fi — > e'j) depends on the combi- 
nation YlY u , which through the Casas-Ibarra parametrization [52] can be expressed in 
terms of an orthogonal matrix R. In the case of hierarchical gauge-singlet Majorana 
neutrinos, one finds that (in the notation of section [2]), 

(Y v )ij « y/MlSoRuiy/K)^. (27) 

Thus, the branching ratio for \x — > e'j depends almost exclusively on one mixing angle 
in R, which we denote as Q\ (we are working under the assumption that R is real). To 
demonstrate the effect of the hidden sector on the branching ratio, we choose a typical 
point in the MSSM parameter space: tan/3 = 20, fi > 0, ao = 0, M = 200 GeV, 
Mi/2 — 760 GeV. In all cases, we set the hidden-sector mass scale at 10 12 GeV. 

We show in Fig. [5] the effect of the hidden sector on the branching ratio for fi — > e'j as a 
function of 9\. In each panel the solid blue curve denotes the prediction for the branching 
ratio in the absence of hidden-sector self-interactions, and the dashed red curve represents 
the prediction with the hidden-sector effect turned on. The top left panel displays the 
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result for 5 7 = 2/3, a 7 = log (5 x 10 11 GeV//x), the top right panel displays the result 
for 6 7 = —2/3, a = log (5 x 10 11 GeV//i), the bottom left panel displays the result for 
6 7 = 2/3, a = log (5 x 10 18 GeV//i), and the bottom left panel displays the result for 
6 7 = -2/3, a = log(5 x 10 18 GeV//i). 




Figure 5: Four panels demonstrating the effect on fi — > e*y for the four behaviours of the 
anomalous dimension outlined in the text. Top left: 5 7 = 2/3, a 7 = log (5 x 10 11 GeV//i), 
top right: 6 7 = —2/3, a = log (5 x 10 11 GeV/fi), bottom left: 6 7 = 2/3, a = log (5 x 
10 18 GeV/fi), and bottom right: 6 7 = -2/3, a = log(5 x 10 18 GeV/fi). 

We see that the largest suppression of the rate occurs for the UV-free anomalous 
dimension displayed in the top two panels. This observation simply reflects that in the 
UV-free setting we chose the pole of the anomalous dimension at 5 x 10 11 GeV, a factor 
of two lower than the hidden sector scale, whilst in the IR-free case we placed the pole 
above the reduced Planck mass. As a result of these choices, the UV-free case has a larger 
anomalous dimension over the integration range. In this UV-free case, the anomalous 
dimension grows during the running of the RGEs from the unification scale to the hidden 
sector scale. Since the seesaw remains active until the Majorana scale is reached, we see 
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that the growth in the hidden-sector anomalous dimension in the UV-free scenario serves 
to influence increasingly the rate for \x — > - in this case suppressing the rate by up to 
an order of magnitude. This observation is consistent with our semi-analytic treatment 
in the previous section. 

We can see how the allowed parameter space changes with the hidden-sector physics 
by examining the allowed parameter regions in the conventional CMSSM M — Mi/ 2 
plane. In Fig. [6] we examine the hierarchical Majorana gauge-singlet neutrino case with 
Q\ = 1.4, displaying the large effect that the hidden sector may have on the predicted 
rate for \i — > ej. The solid white area in each of the three panels indicates the region that 
is currently excluded by the experimental upper limit (BR(/z — ► ej) < 1.2 x 10~ n ). The 
left panel displays the allowed region with the hidden sector effect turned off; the second 
panel shows the effect with the allowed regions for b 7 = 2/3, a 7 = log(5 x 10 11 GeV//x), 
and the third panel displays the effect with 6 7 = —2/3, a 7 = log (5 x 10 11 GeV//i). In 
each case we set the Majorana mass scale to 10 14 GeV and the hidden-sector scale to 
10 12 GeV. 

We see that the hidden sector can change the allowed parameter space quite dra- 
matically. Comparing the last two panels with the first one, we see that some regions 
that were excluded have become allowed and regions that were allowed have become 
excluded. These figures indicate the subtle effects the hidden sector may have. From 
the semi-analytic treatment in the previous section, we saw that for a large enough 
hidden-sector coupling, we expect the rate for ji — > to become suppressed. However, 
as mentioned earlier, the hidden sector also alters the sparticle spectrum, specifically the 
mass eigenvalues of the scalar particles, and this has a nontrivial effect on the branch- 
ing ratio for fi — > ej. In order to determine the full effect of the hidden sector on the 
branching ratio for \x —>■ ej, we require a full numerical treatment as displayed in Fig. [61 
In regions where the hidden sector serves to eliminate parameter space, the change in the 
sparticle spectrum is the dominant effect. 

In Fig. [7J we show the effect of the hidden sector with an IR-free anomalous dimension 
in the same parameter space as Fig. [6j Again, we see that the hidden sector can change 
dramatically the allowed parameter space. The first panel displays the allowed region 
with the hidden-sector effect turned off, the second panel shows the effect on the allowed 
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Figure 6: The allowed parameter space with BR(fi — > ej) < 1.2 x 10~ n in the presence 
of hierarchical gauge-singlet neutrinos and a UV-free hidden-sector theory, with the solid 
white regions denoting excluded regions. Left panel: allowed parameter space without 
hidden-sector dynamics, middle panel: allowed parameter region with & 7 = 2/3, a 7 = 
log (5 x 10 11 GeV//j), and right panel: allowed parameter region with 6 7 = —2/3, a 7 = 
log(5 x 10 11 GeV///). In all cases M R = 10 14 GeV, fj, > 0, tan/3 = 20, a = 0. 

regions for 6 7 = 2/3, a 7 = log (5 x 10 18 GeV//i), and the third panel displays the effect 
with 6 7 = —2/3, a 7 = log (5 x 10 18 GeV///). In all cases we have taken 9i = 1.4 with the 
Majorana mass scale at 10 14 GeV and the hidden sector mass scale at 10 12 GeV. 

4.2 Degenerate Gauge-Singlet Neutrinos 

In the degenerate gauge-singlet neutrino case the combination Y„Yj yields no dependence 
on any particular angle in the orthogonal matrix R. In this case, once Y^Yj is given, 
the branching ratio is fully determined. We begin by displaying the parameter space 
in the M — M1/2 plane. The three panels in Fig. [8] show the changes in the allowed 
parmeter space in the presence of the hidden sector. The first panel displays the allowed 
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Figure 7: Hierarchical gauge-singlet neutrinos with an IR-free hidden-sector theory, dis- 
playing the allowed parameter space with BR({i — > ej) < 1.2 x 10~ n , with the solid 
white regions denoting excluded regions. Left panel: allowed parameter space with- 
out hidden sector dynamics, middle panel: allowed parameter region with 6 7 = 2/3, 
a 7 = log (5 x 10 18 GeV/fi), and right panel: allowed parameter region with 6 7 = —2/3, 
a 7 = log(5 x 10 18 GeV//x). In all cases M R = 10 14 GeV, M hiddcn ~ 10 12 , /i > 0, 
tan/? = 20, and ao = 0. 

region with the hidden sector effect turned off, the second panel shows the effect with 
the allowed regions for b 1 = 2/3, a 7 = log (5 x 10 11 GeV///), and the third panel displays 
the effect with b 1 = —2/3, a 7 = log (5 x 10 11 GeV//i). In all cases we have taken Q\ = 1.4 
with the Majorana scale at 10 14 GeV and the hidden-sector scale at 10 12 GeV. 

In Fig. [9] we repeat the analysis for the IR-free anomalous dimension case with pa- 
rameters 6 7 = 2/3, a 7 = log (5 x 10 18 GeV///). The third panel displays the effect with 
6 7 = —2/3, a 7 = log (5 x 10 18 GeV//t), respectively. Again we see the competition be- 
tween the expected suppression in the rate for \x — > e'j and the alteration of the sparticle 
spectrum. In either case we can see that the hidden sector has a dramatic effect on the 
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Figure 8: The case of degenerate gauge-singlet neutrinos with an UV-free hidden- sector 
theory, displaying the allowed parameter space with BR(fi —>■ ej) < 1.2 x 10 -11 , with 
the solid white regions denoting excluded regions. Left panel: allowed parameter space 
without hidden-sector dynamics, middle panel: allowed parameter region with 5 7 = 2/3, 
a 7 = log(5 x 10 11 GeV/fi), and right panel: allowed parameter region with 6 7 = —2/3, 
a 7 = log(5 x 10 11 GeV///). In all cases M R = 10 14 GeV, M hidden = 10 12 GeV, [i > 0, 
tan (3 = 20, and a = 0. 

rate for [i — > ej and that a knowledge of the hidden sector is not only required to make 
accurate predictions of the low-energy mass spectrum, but is also required to predict the 
level of expected charged-lepton flavour violation. 

5 Hidden-Sector Dynamics and Seesaw Reconstruc- 
tion 

Since the hidden sector affects not only the supersymmetric spectrum but also the 
radiatively-induced seesaw charged LFV, the hidden sector impacts the ability to re- 
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Figure 9: Degenerate gauge-singlet neutrinos with an IR-free hidden-sector theory, dis- 
playing the allowed parameter space with BR({i — > ej) < 1.2 x 10~ n , with the solid 
white regions denoting excluded regions. Left panel: allowed parameter space with- 
out hidden-sector dynamics, middle panel: allowed parameter region with 6 7 = 2/3, 
a 1 = log (5 x 10 18 GeV/fi), and right panel: allowed parameter region with 6 7 = —2/3, 
a 7 = log(5 x 10 18 GeV///). In all cases M R = 10 14 GeV, M hidden = 10 12 GeV, [i > 0, 
tan (3 = 20, and a = 0. 

construct the seesaw from low-energy observations. As we see from eq. fl2Tj) . charged LFV 
fixes the elements of Y^Yj , while low-energy neutrino observations determine the PMNS 
matrix and k in the Casas- Ibarra parametrization [52], eq.(pQ). It has been shown [55] 
that, in principle, in the seesaw extension of the MSSM, a unique map Y^Yj, k — > Y u , M. 
exists, and thus the seesaw parameters can be reconstructed from low-energy data. What 
is important in this reconstruction is that in the seesaw extension of the MSSM all the 
non-seesaw interactions are in principle known from low-energy data, so the RGE mixing 
responsible for the Y„Yj generating charged LFV has only the seesaw parameters as 
unknowns. 
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In our preceding paper [19] , we made the case that, for hidden sectors that could 
be effectively parametrized by a two-parameter characterization of the anomalous di- 
mension of the dominant scalar mass-squared mediation operator (after removing ex- 
ternal line wave-function renormalization effects by rescaling the input parameters at 
M) that the parameters could in principle be fit, and in most cases distinguished from 
other parametrizations by measurements of the low-energy soft supersymmetry-breaking 
terms. This reconstruction can be done using combinations of the soft parameters which 
are unaffected by the seesaw. 

In these circumstances one would still be able to do seesaw reconstruction, as all the 
non-seesaw contributions to the soft-parameter RGEs would be known. However, the 
reconstruction does depend on the assumption that the slepton soft mass matrices at the 
mediation scale M are proportional to the identity, as in the CMSSM case. 

Hence, if the LHC and a linear collider reveal evidence for a distorted superpartner 
spectrum, consistent with strong hidden-sector dynamics, it will be essential to recon- 
struct this dynamics as proposed in [49] first, before one can use other low-energy ob- 
servables, such as charged-lepton flavour violation, to reconstruct the parameters of the 
seesaw mechanism. 

Conversely, if the hidden sector dynamics is sufficiently complicated to defy convenient 
parametrization and experimental determination, then the unknown hidden-sector con- 
tributions to the RGE flow of the soft supersymmetry-breaking parameters will present 
an irreducible obstacle to using low-energy data to determine seesaw parameters. 

6 Comments and Conclusions 

Recent theoretical observations have demonstrated that low-energy predictions in su- 
persymmetric models, such as the supersymmetric spectrum itself, can be significantly 
influenced by a dynamical hidden sector used to break supersymmetry [15J HHJ H7J HH] . In 
this paper, we have examined the effect of a dynamical hidden sector on seesaw induced 
charged-lepton flavour violation. Since both the seesaw sector and the hidden sector sit at 
intermediate scales, an interval exists where both effects are active during the renormal- 
ization group running from the unification scale to the weak scale. The combined effect 
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may alter the usually expected charged-lepton flavour violation in the seesaw extension 
of the MSSM. 

In our analysis, we parameterized the effect of the hidden sector through a simple 
anomalous-dimension Ansatz. In any realistic model, the actual behaviour of the hidden- 
sector coupling would be determined by the theory, but we see that from our simple 
parametrization that we can capture the effect of moderately-coupled IR-free and UV-free 
theories. At a comparison point, we see that the self-coupling in the cubic superpotential 
theory requires A > 5 in order for the hidden sector to have a significant effect on the rate 
for /i — > cy. Our model-independent parametrizations ensured perturbativity in the range 
of integration and, in particular, we ensured the the hidden sector remains perturbative 
up to the reduced Planck scale. In order to generate a large effect on fi — > cy, we need 
a moderately strongly-coupled hidden sector. In the UV-free theory, we placed the pole 
of the anomalous dimension just below the hidden-sector scale itself. This ensured that 
the anomalous dimension becomes large over the range where the seesaw is active. If the 
hidden sector does not become strongly coupled until well outside the interval between 
the hidden-sector and Majorana scales, the effect on flavour violation becomes minimal. 
We can see this effect with the IR-free theory where we place the pole just above the 
reduced Planck mass. Since we started running; the theory at M x = 2 x 10 16 GeV, the 
anomalous dimension started from a smaller value relative to the end-point of the UV-free 
anomalous dimension in our examples. 

Finally, we noted that, in the presence of strongly-coupled hidden sectors whose dy- 
namics distorts significantly the pattern of soft supersymmetry-breaking scalar masses- 
squared, the ability to use soft parameters in a reconstruction of the seesaw depends on 
establishing previously the ability to use the TeV-scale observables to reconstruct the 
hidden-sector dynamics in ways that we have previously analyzed in [49]. Otherwise, the 
observable effects of the hidden sector may impede our view of the seesaw sector. 

Acknowledgements 

We are deeply grateful to Graham Ross for sharing with us the parametrization of the 
hidden-sector renormalization effects presented in Section 4, and also for his encourage- 



20 



ment during the early part of this work. BC and DM would like to acknowledge the 
support of the Natural Sciences and Engineering Research Council of Canada. 



Appendix: One-Loop MSSM Calculation for lj —> / z + 7 

We calculate the rate for \i — > c-f at one loop in the MSSM after running the full sys- 
tem of MSSM RGEs. We follow the notation in [56]; similar formulae can be found 
in |51j . The interactions leading to the LFV process lj — > li + 7 involve the effective 
Lagrangians describing the neutralino-lepton-slepton and the chargino-lepton-sneutrino 
systems. Written in the mass eigenbasis where O ne diagonalizes the neutralino mass 
matrix, Ol and Or diagonalize the chargino mass matrix, XJ~ { diagonalizes the charged- 
slepton mass matrix, and Ug diagonalizes the weak-scale sneutrino mass matrix, we have 

3 4 6 

C = E E E N iabfbEiX Q a + N^/faxl + C. C. (28) 
i=l a=l 6=1 

and 

3 2 3 

C = E E E C i^bE iX - + C^eat, + c c. (29) 

i=l a =l b=l 

where 

NL = -^tanMU^)^ (30) 
NSt = ^ (tan*w (U ? )* (O nc ) al + (U f )* (O nc ) a2 - (U f )^ +3) (O nc ) a3 ) , (31) 



and 



CL = J*^ - (OL) a2 (U fl )*, (32) 
V2mw cos p 

= -ffa (Or) o1 (U fi )i • (33) 

The on-shell amplitude for Zj — > Z, + 7 has the general form 

A* = ee*Zi (p - 5) {im h a^q u (A L L + A R R)) I, (p) , (34) 

where we have used Dirac spinors Z» (p — g) and Z 3 - (p) for the charged leptons z and j with 

momenta p — q and p, respectively; L = (1 — 7 s ) /2 and R = (1 + 7 s ) /2. Each of the 
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dipole coefficients A^ and receives contributions from the neutralino-lepton-slepton 
and chargino-lepton-sneutrino interactions, namely, 



a Jj = aP+a®, 



(35) 



and 



At, - 4 (n) 4- 4 (c) 



(36) 



where A^ , A^' , A^ , A-j^ can be evaluated from the Feynman diagrams in Fig 
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The functions Ji (x), J2 (^), ^3 J a (x) are defined as 

1 — 6x + 3x 2 + 2x 3 — 6x 2 In x 



Ji(x) 
J 2 (x) 
Js(x) 
J 4 (x) 

Finally, the decay rate for Ij 

r ( l 7 ~ 



6(l-x) 4 

1 — x 2 + 2x In x 

2 + 3a; — 6s 2 + x 3 + 6xlnx 

6(l-x) 4 : 
-3 + 4x - x 2 + 21nx 

£j~ + 7 is given by 



C + 7) 



^4(KI 2 + IAn0 



(37) 



(38) 

(39) 
(40) 



(41) 
(42) 
(43) 
(44) 



(45) 



and i — 1, j — 2 for — > e + 7. 
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